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Abstract: The spatial dependence of static meson correlation functions at finite baryon 
density is studied in the Nambu - Jona-Lasinio model. In contrast to the finite temperature 
case, we find that the correlation functions at finite density are not screened but exhibit 
long-range oscillations. The observed phenomenon is analogous to the Friedel oscillations in 
a degenerate electron gas. 



1. Introduction 



Recently, there is a growing interest in studying the properties of hadronic correlation 
functions. The correlation functions form a common ground, where QCD inspired effective 
theories can be tested against each other as well as against lattice results, where available 
i- 

The hadronic correlation functions have been intensively studied in lattice simulations 
of QCD at finite temperatures [0, |^, ^, The results of these simulations show that 
at large distances the correlation functions decrease exponentially. The screening masses 
characterizing this type of behaviour are, at T > T^, degenerate for chiral partners. This 
reflects the restoration of chiral symmetry at high temperatures. Moreover, as it was pointed 
out by Eletskii and loffe @] (see also 0), the values of the screening masses are close to 27rT 
for mesons and SvrT for baryons. Such results are typical for a gas of non-interacting quarks. 
Therefore, the lattice results show that at high temperatures quarks behave, to large extent, 
as free particles. 

The results of the lattice simulations have been interpreted in terms of simple models 
by several authors: Gocksch §1 argues that the hadronic screening masses can be well 
described in the linear sigma model; Hansson and Zahed claim that the static correlation 
functions in high-temperature QCD can be well calculated from an equivalent problem of 
non-relativistic quarks in a dimensionally reduced theory; similarly, V. Koch et al. argue 



that the propagation of a light fermion in a spatial direction at high temperatures can be 
described effectively by a two-dimensional Schrodinger equation. 



In Ref. |TT|| the temperature dependence of the meson screening masses was computed 
and discussed in the framework of the Nambu - Jona-Lasinio (NJL) model. The results are 
in qualitative agreement with the lattice simulations. In the present paper, we explore the 
static meson correlation functions at finite baryon density, using the methods developed in 
1^ |TT]. In this way, we complement and generalize our previous work. We note that, for 
studying systems at finite baryon density one has to rely on models, since lattice results are 
available only for vanishing density. 

The NJL model was introduced already in the early sixties as a theory of interacting 



nucleons [12 . Later it was reformulated in terms of quark degrees of freedom. Numerous 



calculations demonstrate the success of the model in describing hadronic data (for a recent 



review see one of the articles in |13]). Two general characteristics of QCD are described 
by the model, namely the chiral invariance which is spontaneously broken in the ground 
state and the chiral phase transition (the chiral symmetry is restored at high temperature 
or density). On the other hand, the model is not renormalizable and not confining. Thus, it 
is an effective theory, valid only at low energies for phenomena, which are controlled by the 
chiral symmetry of QCD. 

Since the NJL model exhibits a chiral phase transition, it is very instructive to use it 
in order to study the temperature (density) dependence of various physical quantities like. 



e.g., the quark condensate [Q, the pion decay constant |T5|, |1^], and the dynamic masses 
of hadrons |T^, IB, T^. The model has also been used to investigate the thermodynamic 
properties of the quark-meson plasma near the chiral phase transition [|T^. 

In this paper we study the in-medium static meson correlation functions at finite baryon 
density. In contrast to the behaviour known from the high temperature calculations, we find 
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that the correlation functions at finite baryon density are not screened but oscillate in space. 
This behaviour is analogous to the Friedel oscillations in a degenerate electron gas pO[; they 
are caused by the existence of the sharp Fermi surface. 



We find that the NJL approach to the correlation functions breaks down already at 
moderate densities. Thus, to complete the picture we compute the meson correlation function 
at higher densities in perturbative QCD. Keeping only the leading term, we find oscillations 
with a period 5r = tt/pf, where pf is the Fermi momentum of the quark sea. 

The paper is organized as follows. In the next section we define the model. In Section 3 
the details of the calculations are presented. Section 4 contains our results on the density 
dependence of the quark and meson masses. In Section 5 we discuss our results on the 
correlation function. The results based on the perturbative QCD are presented in Section 6. 
We summarize the paper in Section 7. 

2. Definition of the model 



Our calculations are based on the following form of the Lagrangian 



3 Q _ 

a=0 ^ 



Here ■?/' is the Dirac field with additional fiavour {Nf=2) and colour {Nc=3) degrees of 
freedom, aa are the Pauli matrices (with cro=l), Gs is the coupling constant, and m is the 
current quark mass (m^ = = m). 

In the chiral limit, m ^ 0, the Lagrangian (P is invariant under the unitary transfor- 
mations Uy(l), Ua(1) and the chiral symmetry SUi(2)xSUR(2). These are fundamental 
symmetries of QCD, the underlying theory of strong interactions. In the QCD vacuum the 



axial Ua(1) symmetry is broken due to the instanton effects |jT9[ and SUi(2)xSUii;(2) is 
spontaneously broken down to SUy(2). In the Nambu-Jona-Lasinio model the spontaneous 
breaking of the chiral symmetry is reproduced. The explicit breaking of the Ua(1) symme- 
try, accounting for the instanton effects, can be modelled by adding an extra term to the 
Lagrangian ([^). In order to keep things as simple as possible, we do not take this term into 
account here. This does not affect our results, since we do not consider the rj — r]' channel. 

The self-energy of quarks is obtained in the Hartree-Fock approximation by a self- 
consistent solution of the Schwinger-Dyson equation 
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(JflTr [(JaS{p)] - (TaS{p)(Ta 

-Fio-aTsTr [iaa'j5S{p)] - i(yal^S{p)iaal-> 



a=0 ' 



(2) 



Here Tr denotes the trace over flavour, colour and spinor indices, and S{p) is the quark 
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propagator 

S{p)^^ =^ — T, — m + ie. (3) 
Using the last expression we define the quark condensate in the following way 

m = -TjTr S{x = 0- = -AN^i / J" 4 

The sum M = S + m is the constituent quark mass. Eqs. @) and (§) lead to a simple 
relation between M and the condensate 

M = m-2Gs{qq)- (5) 
The zeroth-order correlation function is defined by the expression 

xS(g) = 2^iV,Sp j ^^[V^S{p + Q)VaS{p)]. (6) 

Here A = P {A = S) corresponds to the pseudoscalar (scalar) channel, = {uj,q) is the 
external momentum, Tp = i'j<^, Ts = 1, and Sp denotes the trace over the spinor indices. In 
the random phase approximation the full correlation function has the form 



1 - GsAiQ) 



Xaa{Q) = . r w^^^ - 



It is important to realize that the in-medium correlation function Xaa{Q) depends on 
the variables u;^ and = separately. From now on we shall use the notation xaaI^"^, (f) 
in order to exhibit this dependence explicitly. The dynamic mass is defined by the position 
of the lowest lying pole of XAAiy^"^, 0). In our case it is easily found by solving the equation 



1-G'5X?1«„,0) = 0. (8) 

The screening mass, on the other hand, is defined by the asymptotic behavior of the static 
correlation function in r-space, namely 



. _ ito ^l!^, (9) 



where 



XAAir) = / XAAiO, q') e^^'^ = dqq Xaa{^. q') e^^\ (10) 

— oo 

In Eq. ( p!oD the angular integrals were performed using the fact that Xaa(0, q^) is a function 
of q^ only. In vacuum, the dynamic and screening masses are equal due to Lorentz invariance. 
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On the other hand, at finite temperature or density the two masses can be different, since 
the heat bath introduces a preferable reference frame. The relation between the dynamic 
and the screening masses is discussed more thoroughly in |ll|]. The main purpose of this 
paper is to compute the static correlation functions defined by ([T0|). 

3. Physical quantities at finite density 



i) Imaginary time formalism 



In our calculations we adopt the imaginary time formalism PO] , pT| . Formally, this is 
achieved by the following substitution pT 



(2^ 



/ 27r 
c 



(2vr) 



'IV. 



Here C denotes an integration contour in the complex energy plane The position of the 
contour is fixed by the value of the chemical potential /i, see Fig. 1. In the case when the 
integrand depends additionally on the external frequency, e.g., / = f{p^,p,u), the integrals 
on the right hand side of formula (0) should be evaluated for purely imaginary values of u 
and subsequently analytically continued to real frequencies. 



+100 




Figure 1: Integration contour in the complex energy plane used for the evaluation of the matter parts 



The advantage of the decomposition ( pT]) is that we can separate the so called vacuum 
part from the matter part. The vacuum part of a physical quantity does not explicitly depend 
on the occupation of phase space and reduces at T = = to its vacuum expectation value. 
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On the other hand, the matter part depends exphcitly on the occupation of phase space and 
consequently vanishes in vacuum. 

Using Eqs. (^) and ( p^Ij) we can decompose the condensate into two parts 

{qq) = {qq)vac + {qq)mat- (12) 
In the similar way, using Eqs. @ and (|ll]) we find 



g') = AUQ') + g'). (i3) 



a) Vacuum parts 



In this subsection we recapitulate our calculation of the vacuum parts (more details can 
be found in [TT|). After a Wick rotation of the variables p° = 2^4 and u = iq^, we find 



{qq)vac = -y/l,.ac(M2), (14) 

X%,.ac{-Ql) = h,UM') + qlhMM', ql) (15) 

and 

x1Lac{-ql) = h,vac{M') + {ql + AM')hMM\ ql). (16) 



Here ql = q^ + ql = qi + ■■■ +q\ and the functions Ii^^aci.M'^) and l2,vac{M'^ , ll) defined 
by the expressions below 



h,.ac{M\ql) = -m,j 



(TpE 1 



where d!^pE = d^pdp^. 

Since the integrals (0) and (|T8|) diverge, the quantities (0), (0) and (|l^) are not well 
defined. In order to obtain finite results, we apply a modified Pauli-Villars subtraction 
scheme, where the functions Ii^^adM"^) and l2,vac{M'^ , ll) replaced by 

N 

h,UM') - itaciM') = E Ah,U^l) (19) 

i=0 

and 

TV 

hMM', ql) - itaciM', ql) = ^ Ah,vaMl ql). (20) 

i=0 
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In Eqs. ([19D and (|20|) is the number of subtractions, Aq = 1 and Aq = M. The 
coefficients Ai, for i > 0, have to be chosen in such a way as to provide the finite result for 
Iiyac ^2 vac- Moreover, the correlation functions should be well behaved at infinity (for 
q — » oo), to guarantee the existence of the Fourier transform ([T0|). These requirements lead 
to the following set of constraints 



N 



N 



i=0 



N 



2{N-1) 



0. 



(21) 



Following 0, [TT| one finds 



i=0 



(22) 



and 



27r2 



j=0 




+ ^1 +lnA. 
2A, ' 



(23) 



Substituting Eqs. (p^) and (p3D into (p!4D, ([T5| ) and ([T6|) yields finite expressions for the 
vacuum parts of the condensate and the zeroth-order correlation functions. The correlation 
functions evaluated in such a way are functions of the squared Euclidean momentum. This 
form is convenient for doing the Fourier transform (0). In this case we set g4 = and we 
use Eq. ( pSj) with the substitution q% = q^- 

In order to compute the dynamical masses, we need the function /^^^^(M^, — u;^). This 
can be obtained by performing the analytic continuation of the function defined on the right 
hand side of Eq. ([23|). The substitution ^ iu; ±e (for u; > 0) leads to the following result 
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N ^ 

ltaAM\ -u' ± ^e) = ^ E A <! e(2A, - u) 



W 1 — — arcsm — — 

u V V2AJ V2A, 



+ InAj 



+ e(cj-2A, 



^^(^)-l(a.o.(^).^,.>„A. 



(24) 



in ) Matter parts 

Let us now consider the matter parts of the physical quantities of interest. In the case of 
the condensate we find, using Eqs. (^ and (0), 

{qq)raat = h,raat{M^) , (25) 
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where 



M2' 



(26) 



To evaluate the function Ii^mati^^) collect the contributions from the poles within the 
contour C (see Fig. 1) 



h,mat{M ) = / e{ll - Up) = / - 

TT^ J Ur, TT"^ J 

^ 



u,, 



(27) 



Here 9{x) is the step function, Up — ^p^ + M^, and pp is the Fermi momentum of the 
quarks, i.e., pp — \J [i^ — M'^ . 

As was pointed out in Section 2, the matter parts of the zeroth-order correlation functions 
depend separately on u'^ and g^. Nevertheless, for the extreme cases of a purely timclike 
{Q^ — ("^jO)) ^ purely spacelike momentum {Q^ = (0, q)) the following decompositions, 
analogous to those obtained in vacuum, are possible 



and 



xPp,mati<^', q') = h,mat{M') - h,^,,{M' , u' , q') 



(28) 
(29) 



Here the function l2,mat{M :U , q ) is given by the integral 



l2,mat{M\u'',q^)^AiN,j^j 



dp^ f d^p 



(27r)3 [(p+g/2)2 - M2] [(p-g/2)2 - M2] ' 



(30) 



The calculation of /2,mat(-^^, <^^, ?^) proceeds in the same way as the calculation of the 
function Ii^mat{M'^). Now, as remarked above, the integral over p° is evaluated for imaginary 
u and subsequently analytically continued to real frequencies. The final results (for a; > 0) 
are 



and 



l2,mat{M\u'±ie,0) 



Nc dpp^ 



27r2 
27r2 



/ Up ul-{u±ie)^/A 
log{^-^)-HM^u^±^e) 



(31) 



WM^0„^)=-^7^1n 



2p — q 



2p + q 



2g7r' 



;G{M',q'), 



(32) 
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where ujp — Y^pfT+M^, 



TiM'^.u'^ ±ie) = 



and 



arctan 



= Va;2 _ 4M2 In 



a; 



±i^e{u- 2M) e {2uf - ou) Vou^ - 4M2 



2pF + q 



(a; < 2M) 

{uj > 2M) (33) 



+ 2 

+ Up' 



c^F I In 

2pF - g 

2M2 + gpp + cufV^^ + 4M^ 
2M2 - + urVq' + 4^ 



(34) 



— gin 



M 



We need the analytic structure of the function l2,mat{M'^,0,q^) in the whole complex 
g = |q| plane. To this end we represent l2,mat{M'^: 0, g^) as a sum of two functions, namely 



l2,mat{M', 0, g^) = /i,i(g) + /2y„,(g) 



■(-) 



where 



j(±) 

2,mat 



Pf 



dpp , 2p — 2; ± 26 

m . 

ujp —2p — z ±ie 



(35) 
(36) 



At the end of the calculations we let the infinitesimal e go to zero. The functions l2^at{^) 
have logarithmic cuts parallel to the real axis and stretching from — 2pi7 ± ie to +2pF ± ie. 
At the respective cut the imaginary part of the function l2^at{^) discontinuous. On the 
physical Riemann sheet, this amounts to a change in sign: 



p% + M^-\l^ + M^ 



(37) 



where qr is real and — 2pi,- < < 2pF, 5+ = 2e and 5- — 0. Thus, for the upper sign 
one is above the cut of l2^at ^^<i the lower sign below. Obviously the imaginary part of 

-^2,mat above and below its cut is equal to that of I-z^at above and below its cut. This can 
be summarized by the following equation 



Im/i.mLl^ = gR - ih) = Iml^^^t{z = gR + i5±) 



{+) 



— ImJ. 



(+) f 

2 'mat \ 



(38) 



Consequently, the cuts are arranged in such a way that for z on the real axis, the imaginary 

part of l2tlat{^) cancels that of l2mat{^) and l2,mat{M'^,0,(f) is real for real g as it should 
be. ' 
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Using Eqs. (^) , (^) , (p5|) and (^) we can find the analytic continuation of l2,mat{M'^, 0, g^) 
to purely imaginary values of q. Substituting q = ik we find 

PF 

h„AM\0,-e) = -i^/^^arctau(|). (39) 

' 

The analytic form of l2^mat{M'^,0, —k"^) is easily obtained from (p^ and ( |3^ by analytic 
continuation. 



The integrals (^) and (pO|), which define the functions h^mat and /2,m.at are finite. Con- 
sequently, we are not forced to regularize them. However, if we consider the cutoff to be 
an intrinsic property of the quark-quark interaction, there is no reason not to regularize the 
matter parts in the same way as the vacuum parts. In such an approach, all fermion loops 
vanish when — >■ cxd (or T — > oo) due to a cancellation between the vacuum and matter 
parts, as was discussed in the case of finite temperatures |TT[. This is a desirable property. 



since it guarantees that the condensate and the zeroth-order correlation functions are well 
behaved in the high density/temperature limit. Consequently, we replace the matter parts 
by their regularized counterparts 

N 

h,mat{M') ^ I^^^AM') = E Ah,mat{Al) (40) 

1=0 



and 



N 

l2,mat{M\u\q^) ^ I^^^,,{M\u\q^) = Y.Ah,mat{Aluj',q'). (41) 

i=0 



We note that for /i < 2Ai, which usually includes the region of interest to us, the additional 
terms appearing for i > in Eqs. (|40| ) and (^T]) vanish. Therefore, the regularization of the 
matter parts become effective only at very large densities. 

The structure of Eqs. (0) - (|T6|), (l25|) , ([28|) and (p9D suggests that the following definitions 
will be useful 



and 



/i(M^) = //:,,,(M^) + /^:„,,(M^) (42) 
/2(M^^^g2) = /^_(M^g2 -.;2) + /|;_,(M^u;^g2). (43) 



Eqs. (^) - (|38|), (|39D have a very similar form to the expressions found in [|TT]. The formal 
correspondence can be obtained by the replacement of the functions [exp^Pup) + 1]~^ by 
^9{n — LUp). Nevertheless, there is a crucial difference between the present calculation and 
||1 If : here the cut along the real axis is finite, which leads to a completely different behaviour 
of the correlation functions at large distances. 

4. Density dependence of the quark and meson masses. 

Using Eqs. (|), ([T|), ([l|), (|2|) and (H) we find that 

M = m + MGshiM^). (44) 
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Eq. (44) is the gap equation determining the constituent quark mass M. This equation 
determines M as a function of the chemical potential /i. The dynamic masses of pion and 
sigma are obtained from Eq. (|). Assuming that the gap equation has a non-trivial solution, 
M 7^ 0, we obtain 



m 
M 



+ mdyn,j2{M , m^y^^^, 0) = 



and 



m 
M 



+ {mX^^, - AM'')h{M\ mX^^^, 0) = 0. 



(45) 
(46) 



If the current quark mass is zero, the pion is massless, i.e., a Goldstone boson and the mass 
of sigma is simply 2M. 

We note that Eqs. (B3p and (HBf) are correct not only in vacuum but also at finite density, 
as long as a non-trivial solution to the gap equation exists. Nevertheless, as follows from 
Eqs. ( p4|) and (31), the correlation functions have cuts for arguments larger than 2M. At 



high densities the meson poles merge with the quark-antiquark cuts, which means that in 
this model the mesons can decay into qq pairs. Thus, at these densities there are no isolated 
poles which can be identified with the mass of the pion or the sigma. We circumvent this 
difficulty by defining the mass as the zero of the real part of Eq. (||), i.e.. 



G5RexS(m^^,,0) = 0. 



(47) 



We note that by using Eq. (^) we implicitly neglect the qq widths of the mesons. This 
seems reasonable, since these widths are non-zero only because the Nambu-Jona-Lasinio 
model lacks confinement. 
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Figure 2: Constituent quark mass M as a function of the quark chemical potential 
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In Figs. 2 and 3 we show the density dependence of the constituent quark mass and 



of the dynamic masses of mesons. We use the same values of the parameters as in |11 



Consequently, our results for /i = coincide with the T = results of [|rT|. The values of 
the three regulating masses are: Ai = 680 MeV, A2 = 2.1Ai and A3 = 2.IA2. The coupling 
constant Gs = 0.75 fm^ and the current quark mass m = 8.56 MeV. In vacuum we find the 
constituent quark mass M = Mq = 376 MeV, the pion mass miiyn,TT = 138 MeV, and the 
sigma mass mdyn,a = 760 MeV. The corresponding value of the quark condensate is {qq)^^^ = 
-214 MeV and of the pion decay constant = 94 MeV. The latter can be obtained from 
the Gell-Mann - Oakes - Renner relation: /^m^ = — 2m(gg). 

We note that at T = all values of the chemical potential in the range < /i < Mq are 
equivalent and correspond to the physical vacuum. The reason is that the chemical potential 
fj, is the energy needed to add a quark to the system. Clearly must exceed the minimal 
energy, the vacuum rest mass of a quark Mq, before the Fermi sea of quarks begins to be 
populated. 

In Fig. 2, one can see that the constituent quark mass decreases with the increasing 
value of the chemical potential. At the same time, see Fig. 3, the dynamic masses of pion 
and sigma become degenerate. This fact indicates restoration of the chiral symmetry at high 
density. The overall /i-dependence shown in Figs. 2 and 3 is analogous to the T-dependence 



studied in \11 . 




Figure 3: Pion (solid line) and sigma (dashed line) dynamic masses plotted as functions of the quark 
chemical potential /i. 

The change of the in- medium constituent quark mass M implies, through Eq. (|^), a 
modification of the value of the quark condensate. The calculation of the decrease of the 
condensate with increasing baryon density in the framework of the NJL model is instructive, 
since there are model independent estimates of this quantity, which impose constraints on 
the in-medium behaviour of the condensate. Using the Hellmann-Feynman theorem and the 
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Gell-Mann - Oakes - Renner relation, one finds 



22. 



23] to leading order in the density 



{qq)o mifi 

where Tj^^n is the pion nucleon sigma commutator, (gg)o is the vacuum value of the quark 
condensate, and p is the baryon density. The baryon density is for two flavours given by 
p = 2p|n/37r^, where pf = V f^^ ~ is the quark Fermi momentum. We compute the pion 
nucleon sigma term by using the relation [^] 

-S^AT = = m , (49j 

6 m 

where we assume that the vrA^ sigma term is simply the sum of the irq sigma terms, like in 
the naive quark model. For the parameters obtained above we find T^t^n = 18 MeV, which, 
using Eq. (0), implies that the condensate is reduced by 15 % at the saturation density 
of nuclear matter, po = 0.17 fm~^. In Fig. 4 we show the numerical results for the quark 
condensate (normalized to its vacuum value) as a function of the baryon density (normalized 
to the saturation density). We note that the calculation agrees with the low-density theorem 



(|48|) , when we use the sigma commutator obtained within the model. Since the sigma term 
is much smaller than the empirical value of 45 MeV, the density dependence of the quark 
condensate is too weak. Nevertheless, the fact that Eq. (|48|) is satisfied, shows that the 
calculation is consistent. 




Figure 4: Ratio of the in-medium quark condensate to the vacuum one, plotted as a function of baryon 
density. 
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5. Static correlation functions in space 



i) Analytic structure in complex momentum space 

In this Section we discuss our results on the spatial dependence of static correlation 
functions. By deforming the integration contour, the Fourier transform (|l^) can be rewritten 
as a sum of a few contributions, which are due to the singularities of the correlation function 
in the complex g-plane. The results of Section 3 indicate that Xvaci~Q'^) has cuts (in the 
following we refer to these as the vacuum cuts) along the imaginary axis starting at g = ±2iM 
and going to ±zoo. On the other hand the function Xmat(O) has cuts {matter cuts) which 
are parallel to the real axis. Thus, the full correlation Xaa(— Q'^) function defined by Eq. (|^) 
has at least these two cuts. Moreover, it can have poles for imaginary arguments in between 
the cuts. The general structure of these singularities is shown in Fig. 5. 



upper vacuum cut ^ 


Im q 

2M 

^* 










2pp 









\ 


Req 


matter cuts ^ 


^* 






lower vacuum cut / 









Figure 5: Structure of the singularities of the correlation functions in the complex q — |q| space. 
The position of the pole is determined by 

l-G'5X?i(0,-m^) = 0. (50) 

A pole can appear only in the interval < m* < 2M, since the vacuum cut starts at 
= 2M. In vacuum, due to the Lorentz invariance, we find that m* = rridyn- On the 
other hand, for /i > Mq, ^ T^dyn in general, because the functions l2,mat{M'^,k'^,0) and 
l2,mat{M'^, 0, —k'^) are different; compare Eqs. ( [3T| ) and (pQ]). The contribution of the pole to 
the Fourier transform is of course of the form ~ exp(— m*r). However, there are additional 
contributions from the vacuum as well as the matter cuts. The integral around the vacuum 
cut can be written as follows 
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47r2 



POO p -, 

/ dk k \xaa{0, -k^ + le) - xaa(0, -k^ - le)] e'^^ (51) 



For large values of r, the exponential factor in (^) cuts down the integrand very quickly as 
a function of k. Thus, the remaining factors, which are slowly varying, can be approximated 
by their value at ~ 2M. Consequently, for r — ^ oo the contribution of the vacuum cut is 
also of the exponential form, ~ exp(— 2Mr). The contribution of the matter cut is 

+2pF 

XAAir) = ^ J dqq [xaa(0, [q + tS.]') - x^a(0, \q + e^''^ (52) 

-2PF 

where 6± is defined below Eq. (^7|) . The imaginary shift of the argument, denotes 
that the function should be calculated just above/below the matter cut. We note that in 
contrast to the finite temperature calculation |jll|, the contribution by the matter cut is 
in the present case represented by an integral over the finite range in q. This leads to a 
qualitatively different behaviour of the Fourier transform. 

n) Results 

Let us start the presentation of the results by discussing the zero density case (0 < yU < 
Mo). This is equivalent to the zero temperature case (T = 0) studied in []1T|. In vacuum 



the matter cut vanishes and we are left with only one or two singularities which contribute 
to the Fourier transform. This allows for the simple analysis of the asymptotic behaviour of 
the correlation functions. 

In the case of the pseudoscalar channel we find an isolated pole and a cut. Because of 
the Lorentz invariance, the pole coincides with the dynamical mass. For large values of r 
the contribution from the pole, ~ exp(— m^y^ ^r^); dominates over that of the vacuum cut, 
~ exp(— 2Mor). Consequently, the static correlation function in the pseudoscalar channel 
exhibits a pure exponential behaviour, whose decay rate is characterized by the pion screening 
mass: mscr,n = fndyn,n = 138 MeV. 

In the case of the scalar channel we do not have an isolated pole and the only contribution 
to the Fourier transform is due to the vacuum cut. The latter also leads to the exponential 
behaviour. In consequence the sigma screening mass equals 2Mo and is a little bit smaller 
than the dynamic one computed using Eq. (^T]): mscr,a = 752 MeV and mdyn,a = 760 MeV. 

Our numerical calculations, done in the interval 0.5 fm < r < 6 fm, confirm the re- 
sults of the analysis carried out above: the correlation functions in both pseudoscalar and 
scalar channel decays exponentially with distance, although the numerically evaluated sigma 
screening mass turns out to be slightly larger than 2Mq. The discrepancy between the nu- 
merical and analytical results is due to the fact that numerically we are not able to compute 
at truly asymptotic distances. 

Let us now discuss the finite density calculations. In the case of the pseudoscalar channel, 
we did the numerical calculations in the energy range Mq = 376 MeV < /i < 415 MeV. In 
Fig. 6 we show the corresponding correlation for fi = 376,400 and 410 MeV. In contrast to 
the exponential decay found in vacuum (solid line), the correlation function at finite density 
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(dashed and dotted lines) oscillates, with a power-law decay of the amplitude. Its period 
decreases with the increasing value of the chemical potential. From the physical point of 
view, such oscillations are caused by the existence of a sharp Fermi surface. 
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Figure 6: Correlation function in the pseudoscalar channel plotted for three different values of the chemical 
potential: = 376 MeV (solid line), fi = 400 MeV (dashed Hne), and fi = 410 MeV (dotted hne). 

The oscillatory behaviour of the correlation functions at finite density is well known for the 
non-relativistic degenerate electron gas, where the phenomenon is called Friedel oscillations 



20| . These long range oscillations lead to many interesting phenomena like, e.g, broadening 
of the nuclear magnetic resonance lines. A characteristic feature of the Friedel oscillations (at 
very large distances) is their period Sr = ir/pp- In our case the situation is similar; for very 
large values of r, the contribution from the matter cut dominates over the other ones. Since 
the cut extends over a finite range in q {\q\ < 2pp) the correlation function oscillates with 
the period Tc/pp (at r oo). Nevertheless, the numerical study of the correlation functions 
at very large distances is difficult because the amplitudes decrease with r. Consequently, we 
restrict the numerical calculations to the interval 0.5 fm < r < 6 fm. 

So far we have discussed the results for fi < 415 MeV. For larger values of the chemical 
potential (corresponding roughly to p > |po) the pseudoscalar correlation function acquires 
additional singularity on the real axis for q K,2pp. In the scalar channel, a similar singularity 
appears at even smaller densities (p > |po)- These singularities indicate instabilities of the 
ground state, which lead the system to states of lower energy. However, most likely the 
singularities are artifacts due to regularization procedure and do not correspond to physical 
instabilities. In order to check this point we have redone the calculation using a different 
regularization scheme with a three-dimensional cutoff. Using the parameters obtained in 
Ref. |T^, we find no singularities on the real axis. This shows that these singularities are 



unphysical, since their presence depends on details in the formulation of the model. Thus, 
our regularization scheme can be used only at low densities. 
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On the other hand, the three-dimensional cutoff regularization is unsatisfactory because 
it exphcitly breaks Lorentz invariance. This obscures the relation between the correlation 
functions in timelike and spacelike regions. Moreover, in the three-dimensional cutoff scheme, 
one does not recover the well known screening of the correlation function at finite tempera- 
tures. Instead, the correlation function oscillates, much like the Friedel oscillations at T = 0, 
with a period vr/A, where A is the momentum cutoff. We stress that this behaviour is an 
artifact, due to the finite range of momenta (|g| < A) available in the Fourier transform. 

Consequently, we stick to the covariant Pauli-Villars method, in spite of its shortcomings. 
However, since the results are inconclusive at high densities, p ~ Po and higher, we restrict 
the calculations to small densities p < |po- The physics at very large densities will be 
addressed in the following section using a different (although closely related) approach. 

6. Perturbative QCD calculation 



The formalism developed so far can also be used to study correlation functions in per- 
turbative QCD. In this case the leading term in the meson correlation function is the lowest 
order quark loop contribution. Thus, at large densities the correlation function in r-space is 
given by 



1 °9 

X?i(r) = ^ j dqqxfAO,q')e''''. (53) 

— oo 

The Fourier transforms of the vacuum XAA,mc(0; 0'^) matter XAA.mati^jQ"^) parts define 
the functions XAA,vacir) and XAA,matir), respectively, with 



X?i(r) = xiVcW + xSl„«.W- (54) 

The Fourier transform (|53|) can be computed analytically. Using the methods developed 
in [^, we find that 

N ^ 

A,.acir) = ^ E AA? m,i2A,r) + 2K,rK,{2K,r)] . (55) 



where Ki and K2 are the modified Bessel functions pj]. The matter piece, on the other 
hand, is given by 



xPp,.natir) = f 4 - 4^ + [Gi(r) + G,(r)] , 

47T f \ 1" ly ^ (J If If (J If j 



(56) 



where 



^ 2 io ^/^A2 ' " ■ S 

and 



= ff ^J§=|-(2p.) = -i:^A.AM2A.r) (57) 



d Ai f°° dp , , 

G2{r) = --Y.^ -^=^cos(2pr). (5^ 

or I Jpp Jp2 J^2 
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The properties of the modified Bessel functions ||25] imply that the vacuum part (|55|) is 



exactly canceled by the Gi{r) term of the matter part (|56D. The cancellation is quite general 
and in fact independent of the number of subtractions N and of the values of the regulating 
masses Aj. 

The vacuum part of the correlation function Xaa vad^^ ^f) is divergent and must be regu- 
larized. Nevertheless, the Fourier transform x'aa vacix) remains finite for r > also when one 
sends the cutoff masses to infinity. To see how this works we consider the Fourier integral 
of the regularized vacuum correlation function. This is well defined, so that the integration 
contour can be deformed into a path around the cut on the positive imaginary axis starting 
at p = 2iM, where M is the quark mass. The integral around the cut converges for r > and 
the contributions of the regularization terms, which correspond to cuts starting at p = 2iAj 
are exponentially suppressed. Consequently, the role of the regularization terms is to make 
the Fourier integral finite and thus to make the deformation of the contour possible. After 
the contour has been deformed, the integral is finite, for r > 0, even when one sends the 
cutoff masses to infinity. In the limit Aj ^ oo and M — > 0, we find 

G,(r) = ^^^^ (59) 

The form of G2{r) implies that the correlation function Xpp(^) is not screened but oscillates 
in space with a period 6r = tt/pf- For massless quarks, the behaviour of X^ssi^) is identical, 
since the two channels are degenerate in this limit. 

7. Summary 

In this paper, we have studied the structure of static meson correlation functions at finite 
baryon density within the NJL model. We have restricted our work to the pseudoscalar and 
scalar channels. 

In general the vacuum correlation functions in the meson channels are screened and, due 
to Lorentz invariance, the screening mass of a stable meson equals its dynamical mass. For 
unstable mesons the screening mass equals the mass of the branch point, since the asymptotic 
form of the correlation function picks out the lowest lying singularity. As demonstrated in our 



previous work [11], these results are qualitatively reproduced by the Nambu-Jona-Lasinio 
model: the pion screening mass is equal to its dynamic mass, whereas the sigma screening 
mass is 2Mq, where Mq is the constituent quark mass in vacuum. [The latter result is 
of course an artifact of the NJL model, which lacks confinement and of the approximation, 
which lacks a coupling of the sigma to the two-pion continuum. In a more realistic treatment 
the screening mass in the sigma channel should equal 2171.,^.] At finite temperature and zero 
density, the correlation functions are again exponentially damped, although the screening and 
dynamic masses differ. However, as we have demonstrated here, the correlation functions at 
T = and finite density differ qualitatively from those in vacuum: they exhibit long ranged 
oscillations, of the Friedel type, rather than exponential damping. 

In order to understand this effect, we carefully explored the analytic structure of the 
correlation functions. We found that the appearance of the oscillations is connected with the 
existence of a cut of finite range in the complex momentum plane. This cut is responsible for 
the leading contribution to the correlation function at large distances. The length of the cut. 
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which is proportional to the Fermi momentum of the constituent quarks, is reflected in the 
oscillation period at large distances 6r = ir/pp. Consequently this form of the correlation 
function is quite general and is expected in all normal Fermi liquids. In particular, the 
existence of the oscillations is independent of whether the basic fermionic degrees of freedom 
are quarks or nucleons. Therefore, we feel the oscillatory behaviour of the correlation function 
at finite density will not change qualitatively by confinement. 

The fact that the correlation function exhibits Friedel type oscillations indicates that 
at finite density it is impossible to obtain the information about timelike excitations by 
studying spacelike correlation functions. Their long distance behaviour is dominated by low 
lying particle-hole excitations, and consequently not connected with the dynamical mass. 



We find, in agreement with the low density theorem p2| , pB| , that the quark condensate is 



reduced as the density is increased and eventually almost vanishes, implying that the chiral 
symmetry is restored at high densities. This, and the fact that the dynamical masses of 
chiral partners (here vr and a) become degenerate with increasing density, seems to indicate 
that the NJL model offers a consistent description of the restoration of chiral symmetry 
at finite baryon density. Nevertheless, the spurious singularities we encountered in the 
correlation functions, suggest that the model in the present formulation cannot be naively 
used for densities p > |po- Therefore, we complemented our investigation with a calculation 
in perturbative QCD, which again indicates the presence of oscillations also at very high 
densities. 
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